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ABSTRACT 

We present multi- wavelength imaging and near-IR spectroscopy for ten gravitationally lensed galaxies 
at 0.9 < z < 2.5 selected from a new, large sample of strong lens systems in the Sloan Digital Sky 
Survey (SDSS) DR7. We derive stellar masses from the rest-frame UV to near-IR spectral energy 
distributions, star formation rates (SFR) from the dust-corrected Ha flux, and metallicities from the 
[N II] /Ha flux ratio. We combine the lensed galaxies with a sample of sixty star- forming galaxies from 
the literature in the same redshift range for which measurements of [N II] /Ha have been published. 
Due to the lensing magnification, the lensed galaxies probe intrinsic stellar masses that are on average 
a factor of 11 lower than have been studied so far at these redshifts. They have specific star formation 
rates that are an order of magnitude higher than seen for main-sequence star-forming galaxies at 
z ^ 2. We measure an evolution of 0.16 ± 0.06 dex in the mass-metallicity relation between z ~ 1.4 
and z ~ 2.2. In contrast to previous claims, the redshift evolution is smaller at low stellar masses. We 
do not see a correlation between metallicity and SFR at fixed stellar mass. The combined sample is in 
general agreement with the local fundamental relation between metallicity, stellar mass and SFR from 
Mannucci et al. (2010, 2011). Using the Kennicutt-Schmidt law to infer gas fractions, we investigate 
the importance of gas inflows and outflo ws on the shape of the mass-metallicity relation using simple 
analytical models. This suggests that the Maiolino et al.| (2008l calibration of the [N II] /Ha flux ratio 
is biased high. We conclude that both an absolute metallicity calibration and direct measurements of 
the gas mass are needed to use the observed mass-metallicity relation to gain insight into the impact 
of gas flows on the chemical evolution of galaxies. 

Subject headings: galaxies: high- redshift, galaxies: evolution, strong gravitational lensing 



1. INTRODUCTION 

Our current understanding of the build-up of a galaxy's 
stellar mass and metallicity over cosmic time involves 
a subtle interplay of gas accretion, star formation and 
galactic winds. The accretion of metal-poor gas from 
the intergalactic medium (IGM) fuels the star forma- 
tion; subsequent generations of stars contribute metals 
to the surrounding interstellar medium (ISM) and power 
galactic-scale outflows of enriched gas. This simple pic- 
ture of chemical evolution is generally accepted, but the 
detailed physical mechanisms behind each of its compo- 
nents are far from well understood. There is as of yet no 
consensus on how to incorporate physically accurate pre- 
scriptions of gas inflows and outflows into semi-analytical 
and numerical models of galaxy formation. 

Observations have not been able to provide much guid- 
ance. In recent years, spectral line measurements have 



at least comparable to and possibly many times greater 
than the star formation rate for the highly magnified 



& Lewis 


1996 


escape t 


le gra 



al. 2003 Weiner et al. 2009 



at all redshifts (Shapley et 


Steidel et al. 


2010). How- 



ever, the extent to which these outflows attect the chem- 
ical evolution of galaxies depends on the mass outflow 
rate, which has only been measured for one high red- 
shift source: Pettini et al. (2002) found a mass loss rate 



1996). It is unclear whether these outflows will 
le gravitational potential of their host galaxy or 
rain back down. The high metallicity of the diffuse IGM 
in galaxy groups and clusters is indirect evidence that at 
least some fraction of the metals in these winds is lost to 
the environment (e.g. Renzini et al. 1993). 

The accretion of metal-poor gas is generally thought to 
happen along filaments and is refered to as cold-stream 
accretion or cold flows. It has received much recent at- 
tention in numerical simulations of galaxy formation as 
the dominant mechanism which fuels the star formation 
of high-redshift galaxies (e.g. Dekel et al. 2009). The ex- 
istence of massive stellar populations at 2; ~ 2 with old 
ages, high star formation rates (SFR), and continuous 
or slowly decreasing star formation histories requires a 
replenishment of their gas reservoirs by accre tion to sus- 
tain this vigorous, exten ded star formation fShapley et 
al. 120051 IPapovichet al.|2006; Erb 2008). However, direct 
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observations of cold flows are still lacking. Line profile 
measurements at z 2 — 3 have found little evidence 
of infalling material, though the expected narrow veloc- 
ity range of infiowing gas can easily be disguised by the 
several hundred km /s spread observed in most outflows 
( Steidel et al.||2010[ ). Additionally, the covering factor of 
filaments for each galaxy is small, such that projection ef- 
fects strong ly reduce the proba bility of detecting them in 
absorption (Steidel et al. 2010). The recent discovery of 
two extremely metal-poor Lyman limit systems at z ^ 3 
in the spectra of bright quasars presents the flrst obser- 
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vational hints of cold flows at high redshift ( Fumagahi et 
aI] |20Tl1 ). 



Until inflows and outflows can be observed directly and 
their mass flow rates can be quantifled, measurements of 
gas-phase metallicities and their relation to stellar mass 
can provide indirect insights into the impact of these gas 
flows on the chemical evolution of galaxies over time. In 
the closed-box model of chemical enrichment, the metal- 
licity uniquely depends on the stellar mass as a record 
of the accumulated conversion of gas into stars. Devi- 
ations from this closed-box model prediction point to- 
wards metal exchanges betw een the galaxy and its envi- 



ronment through gas flows. Lequeux et al. ( 1979 1 first 



noted the increase in metallicity for more massive galax- 
ies. The correlation w as quantified for the S loan Digital 
Sky Survey (SDSS) by [Tremonti etall ( |2004| , who found 
a tight (±0.1 dex) mass-metaliicity relation (MZR) with 
a relatively steep slope, which flattens above 10^°'^ M©. 
A study of dwarf galaxies has extended the relation to 
much lo wer stellar mass es (lO^'^ ^q) with a similar low 

A clear evolution of the MZR 



scatter ( |Lee et ah 2006[ ) 
with redshift has been 



with redshift has been found using galaxy sam ples at 



0.4 < z < 1.0 ( [Savaglio ct al. 2005]), 1.0 < z < 1.5 ([Shap^ 



(Maiolinoet al.||2008 Mannucci et al. 



PI , ana z 

3DD9 ) : at a flxed stellar mass the nietallicity decreases at 
earlier times. 

Much effort is being invested in using both semi- 
analytic models and cosmological hydrodynamic sim- 
ulations to explain the observed MZR. One common 
idea states that galactic winds can more efficiently re- 
move metals from the shallower potential wells of low- 
mass galaxies, thus reproducing the observed MZR (e.g. 
Kobayashi et al. 2007; Scannapieco et al. 2008, Spi- 
toni et al. 2010). In an alternative scenario, low-mass 
galaxies convert their gas reservoirs into stars over longer 
timescales than more massive galaxies. They subse- 
quently show larger gas-to-stellar mass fractions and are 
simply less enriched. Some studies have required galactic 
winds to regulate this mass-dependent star formation ef- 
ficiency (e.g Brooks et al. 200 7, Finlator & Dave 2008). 
In contrast, Tassis et al. (2008) employ a critical density 
treshold for the onset of star formation and consequently 
do not require outflows to reduce the star formation ef- 
ficiency of low-mass galaxies. Finally, a net dilution of 
the gas-phase metallicity through inflows of metal-poor 
gas when the timescale for star formation falls below the 
accretion timescale can explain the observed MZR (e.g. 
Dalcanton et al. 2004). 

Compared to this extensive theoretical effort, obser- 
vational studies of the MZR are much less developed. 
Beyond z > 1, metallicity measurements of individual 
galaxies are still fairly rare. A large investment of near- 
IR spectroscopy is required to detect the rest-frame opti- 
cal nebular emission lines used as metallicity indicators. 
Faint, low-mass galaxies are even more expensive to ob- 
serve, limiting current z ^ 1 — 2 studies to stellar masses 
> 10^-^ M0. Gravitational lensing provides unique op- 
portunities to extend metallicity studies to lower stellar 
masses. Strongly lensed galaxies are typically magnified 
by factors of a few up to > 50 and probe systems that 
are intrinsically an order of magnitude less massive with 
significantly increased observing efficiency. 



Only a handful of the known lensed galaxies in the lit- 
erature hav e robust measurements of stellar mass and 
meta llicity (Teplitz et a D'2000', 'Lemoine-Bu sserolle et 

Bian ct al. 201^. Many 



Mainline et al.j 2009 



,al. „2003[ 

lensed galaxies are discovered serendipitously, and often 
not much follo w-up is invested beyo nd the discovery pa- 
per. Recently Richard et al. ( 2011[ ) presented metallic- 
ity measurements for ~15 gravitationally lensed sources 
from the Hubble Space Telescope (HST) archive. To cre- 
ate a large sample of highly magnified lensed galaxies, 
we have systematically surveyed the galaxy clusters in 
the Sloan Digital Sky Survey (SDSS) Data Release 7 
and identified a uniformly selected, complete sample of 
^150 confirmed lens systems (Bayliss et al. 2011, M. D. 
Gladders et al. 2012, in preparation). With an aver- 
age g-band magnitude of m^B=22.0, the lensed sources 
are sufficiently bright for multi-wavelength photometric 
and spectroscopic follow-up. The stellar populations and 
physical conditions of two galaxies from this sample at 
z — 2.76 and z — 2.92 have been studied extensively with 
multi-wavelength ph otometry and near-IR spectroscopy 
(Wuyts et al. 2012). In this paper we present simi- 
lar observations tor an additional ten lensed galaxies at 
0.9 < z < 2.5. We compare our results to measure- 
ments from the literature to extend the MZR for indi- 
vidual galaxies at z = 1 — 2 down to lower stellar masses. 

The paper is organized as follows. f|2] describes our 
sample selection, observations and data reduction. We 
derive the galaxy properties - stellar mass, star forma- 
tion rate, and metallicity - in S|3] Q summarizes the 
additional data taken from the literature and f|5] presents 
the combined results. We discuss observational con- 
straints on the origin of the MZR in ^ Throughout 
the paper we adopt a fiat cosmology with — 0.3 and 
Hp = 70kms~ ^ Mpc~^. We assume a Chabrier initial 
mass function ( Chabrier|2003 ) and quote all magnitudes 
in the AB system. 

2. SAMPLE SELECTION, OBSERVATIONS AND DATA 
REDUCTION 

The lensed galaxies discussed in this paper are drawn 
from the subset of 26 SDSS lens systems that have 
been targeted by the Gemini Multi-Object Spectrograph 
(GMOS) to measure velocity d ispersions and dyn amical 
masses for the lensing clusters ( Bayliss et al.|20lT ). This 
effort has identified 69 unique lensed background sources 
with redshifts as high as z = 5.2. We obtained near- 
IR spectra for a subset of 10 lensed galaxies within 8 
of the lens systems, mostly targeting the primary giant 
arc in each system. These sources were selected from 
the larger sample based on observability, rest-frame UV 
surface brightness and the size of the arc, so as to take 
optimal advantage of the 42" long Keck/NIRSPEC slit. 
Throughout t he paper we adopt the naming convention 



introduced in Koester et al] ( |2010 ), which designates 

IMM55+DDMMSS. 



lensed sources as SGAS JHH 



2.1. 



Imaging 

The optical imaging of our sample consists of gri 
GMOS pre-imaging obtained to facilitate mask design 
for th e GMOS spectroscopy presented in [Bayliss et al. 
(2011 1. Each image is a stack of two dithered exposures 
of 150 s each, taken with the detector binned 2 x 2 for 
a scale of 0.1454" pixel^^ and reduced using the Gem- 
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Table 1. Photometry 



source 


9 


r 


i 


z 


J 


H 


3.6 ^im 


4.5 /im 



SGAS 
SGAS 
SGAS 
SGAS 
SGAS 
SGAS 
SGAS 
SGAS 
SGAS 
SGAS 



J085137+333114 
J091538+382658 
J103842+484926 
J103844+484914 
J111124+140901 
J113809+275445 
J113809+275439 
J120924+264052 
J134334+415509 
J142039+395454 



21.66±0.06 
22.21±0.08 
21.26±0.07 
21.27±0.11 
21.75±0.11 
21.75±0.07 
21.31±0.11 
20.96±0.06 
20.46±0.08 
21.11±0.07 



21.40±0.05 
21.95±0.07 
20.66±0.06 
20.74±0.11 
21.57±0.11 
21.24±0.07 
21.01±0.11 
20.67±0.06 
20.27±0.06 
20.63±0.06 



21.04±0.05 
21.52±0.07 
20.37±0.06 
20.02±0.11 
21.87±0.11 
21.05±0.07 
20.62±0.11 
20.02±0.06 
20.39±0.07 
20.27±0.06 



21.44±0.09 20.71±0.30 



17.64±0.12 
18.67±0.31 



19.64±0.17 



20.30±0.07 19.96±0.14 



19.13±0.12 
18.68±0.17 
16.83±0.15 
17.52±0.21 
21.12±0.21 
20.19±0.16 
20.52±0.24 
19.32±0.15 
19.32±0.15 
18.75±0.23 



18.87±0.12 

16.16±0.15 
17.08±0.21 
21.06±0.21 
19.95±0.15 
20.71±0.23 
19.69±0.15 

18.61±0.23 



Note. — Magnitudes are quoted in the AB system. 



ini IRApi package. SGAS J091538+382658 and SGAS 
J134334+415509 have additional 300 s GMOS z-band 
data, obtained for the study of two lensed Lyman-g emit- 
ters at z ^ 5 discovered in these fields (Bayliss et al. 



2010|. 

We have carried out near-IR imaging for three lensed 
galaxies with the Near-Infrared Camera and Fabry- 
Perot Spectrometer (NIC-FPS) on the 3.5 m tele- 
scope at the Apache Point Observatory (APO). SGAS 
J120924-K264052 was observed for 3150 s in the H-hand 
on 2008 May 19 UT, SGAS J091538-I-382658 was ob- 
served for 7200 s in the J-band on 2009 February 6 
UT and SGAS J134334+415509 was observed for 5400 s 
in the J-band on 2009, March 16 UT. All data was 
taken under sub-arcsecond seeing conditions. A custom 
pipeline of standard IRAF tasks was used to sky-subtract 
and stack the dithered NIC-FPS images. 

Warm Spitzer observations with IRAC at 3.6 and 
4.5 /im were obtained for a larger sample of ^ 100 
lens systems through Spitzer program 70154 (PI: M. D. 
Gladders). The lensed galaxies studied in this paper 
all have an integration time of 1200 s in each band, 
with the exception of SGAS J091538-f 382658 and SGAS 
J134334+415509. For these two sources, warm IRAC ob- 
servations at 3.6 fim were obtained through Spitzer pro- 
gram 60158 (PI: M. D. Gladders) with a tot al integ ra- 
tio n tim e of 1800 s and 2700 s respectively (Bayliss ct 
al l2010[ ). AU IRAC data was reduced with the MOPEX 
software distributed by the Spitzer Science Center and 
drizzled to a finer pixel scale of 0'.'6pix~^. 

Photometry of the lensed galaxies is performe d using 
largely the sa me method as described in detail in |Wuyts] 



et al. (2010). In summary, we define ridgelines along 
the arcs and create non-circular photometric apertures 
as isophotes of the convolution of these ridgelines with 
the empirical PSF of the images. Final magnitudes are 
measured at an equivalent radius of twice the FWHM 
of the image and aperture-corrected to an equivalent ra- 
dius of 6" based on the curve of growth of PSF refer- 
ence stars. Aperture corrections for the IRAC data are 
applied as detailed in the IRAC handbook^; a different 
correction is made for point sources (< 8") and extended 

^ IRAF (Image Reduction and Analysis Facility) is distributed 
by the National Optical Astronomy Observatories, which are op- 
erated by AURA, Inc., under cooperative agreement with the Na- 
tional Science Foundation 



sourc es (> 8"). We use GALFIT version 3.0 ( |Peng et aL 
20101 to mask neighboring galaxies that fall within the 
photometric apertures and nearby bright cluster galaxies 
with haloes that extend out to the location of the arc. 
This is especially important for the IRAC data due to 
the extended, non-circular PSF of the instrument. Final 



magnitudes are corrected for galactic extinction ( Schlegel 
et al. 1998 ) and presented in Table [T] Photometric un- 
certainties include Poisson noise, absolute zeropoint un- 
certainties, uncertainties from the Spitzer aperture cor- 
rections as detailed in the instrument handbooks, and 
uncertainties from the masking of neighboring galaxies. 
This last effect often dominates the total uncertainty. 

2.2. Near-IR spectroscopy 

We obtained near-IR sp ectra for this sample u sing the 
NIRSPEC spectrograph ([ McLean et ~ar||1998[) on the 
Keck II telescope on the following nights of 2011: March 
19 UT (full night), March 22 UT (hah night), and March 
23 UT (half night). Cirrus was present on the second 
and third nights. The seeing was measured at telescope 
re-focusing, and monitored from the CFHT seeing mon- 
itor; the measured seeing was highly variable, ranging 
from 0.7" to 2.4". We used NIRSPEC's low-resolution 
mode with the 0.76" x 42" slit, and set the position an- 
gle for each source to maximize coverage of the extended 
arc. Figure [T] shows the location of the NIRSPEC slit on 
composite gri images of the lensed gal axies (created from 
GMOS pre-imaging discussed in ^2.1). Each source was 
acquired by offsetting from the brightest cluster galaxy 
on the near-IR slit-viewing camera; source acquisition 
was verified by direct imaging with this camera. During 
spectroscopic integration, the observations were nodded 
along the slit in an ABBA pattern, with typical exposure 
times of 600 s. To enable telluric calibration and flux- 
ing, immediately before or after each target an AOV star 
was observed, ideally located less than 10 degrees from 
the target and with a similar airmass {S{secz) < 0.1). 
Table [2] summarizes the spectroscopic observations. 

The spectra were reduced with the nirspec-reduce 
package written by G. D. Becker. Th e reduction is de- 



scribed in more detail in ' Rigby et al.| \2Qll) and can 
be summarized as using lamp exposures to flatten the 
data, using skylines to calibrate wavelengths, and per- 

^ http:/ /ssc. spitzer. caltech.edu/irac/iracinstrumcnthandbook/ 
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Fig. 1. — Color-composite images of our lensed galaxies from the gri CMOS pre-imaging, 50" on a side. North is up and East is left. 
The 0.76" X 42" NIRSPEC slit is overplotted for each source. For simplicity we show the slit centered on the source; in fact the source 
was placed on the left half and then the right half of the slit for the ABBA nod pattern, with the nods separated by 10-15". 



Table 2. Spectroscopic Observation Log 



source 


PA 


filter 


grating 


A 


t 


AOV star 


extraction 








angle 


[/im] 


[s] 




method 


SGAS J085137-f 333114 


357 


N-6 


33.50 


1.60-2.02 


1800 


HD71906 


each exposure 


SGAS J091538-I-382658 


5 


N-5 


36.72 


1.51-1.79 


5400 


HD71906 


each exposure 


SGAS J103842-I-484926 


49 


N-6 


35.00 


1.93-2.32 


2400 


HD99966 


each exposure 


SGAS J103844-I-484914 


40 


N-3 


34.08 


1.16-1.37 


3600 


HD99966 


stacked nods 


SGAS J111124-I-140901 


77 


N-6 


34.50 


1.81-2.23 


3000 


69Leo 


each exposure 


SGAS J113809-I-275445 


87 


N-5 


36.40 


1.47-1.74 


6000 


HD105388 


stacked nods 


SGAS J113809-f 275439 


87 


N-3 


34.08 


1.17-1.38 


2400 


HD105388 


each exposure 


SGAS J120924-I-264052 


163 


N-4 


35.23 


1.30-1.59 


3400 


HD105388 


each exposure 


SGAS J134334-I-415509 


146 


N-6 


34.50 


1.81-2.23 


2400 


HD109615 


each exposure 


SGAS J142039-f 395454 


312 


N-6 


34.50 


1.80-2.22 


2400 


HD116405 


each exposure 



Note. — Columns are source name, position angle of the slit, NIRSPEC filter name, cross-disperser angle, wavelength coverage, total 
integration time, AOV star used for telluric correction and extraction method. 



forming sky subtraction following Kelson (2003). We 



employ two different extraction strategies, depending on 
target brightness. For an object with fairly bright emis- 
sion lines, we measure the spatial profile from Ha on each 
exposure and use this to optimally extract the spectrum 
from that exposure. For the faintest objects, we stack 
the 2D sky-subtracted exposures for the A and B nod 
positions, fit the spatial profile of Ha and optimally ex- 
tract the stacked spectrum for each nod position. For 
each object, the extraction method is listed in Table [2j 
When the emission is concentrated in individual "knots" 
within the slit, we separately extract and analyze the 
spectrum of each knot. 

Each extracted spectrum is corrected for telluric ab- 
sorption and fluxed based on a AOV standard star spec- 
trum, using the tool xtellcor_general (Vacca et al. 2003). 
Fluxing is thus appropriate for the fraction of the galaxy 
falling inside the slit, not for the whole galaxy. Given 
the highly variable seeing, the absolute fluxing may have 
significant systematic errors. For objects with a contin- 



uous arc morphology, we combine the individual flux- 
calibrated spectra from each exposure with a weighted 
average to produce a single fluxed spectrum, along with 
its 1 a error spectrum. When objects are composed of 
discrete emission "knots" or multiple images, we create 
a single fluxed spectrum for each knot or image in the 
same manner. 

We derive fluxes for the Ha, [N II] and [S II] 
emission lines as follows using the IDL Levenberg- 
Marq uardt least-squares code MPFITFUN ( [Markwardt 
20091. First, to determine the redshift, we ht a single 
Gaussian to the Ha emission; the uncertainty in the mea- 
sured wavelength centroid is used in subsequent fitting. 
Second, we simultaneously flt three Gaussians to the sys- 
tem of the [N II] emission lines and Ha, with the following 
physically motivated constraints: 

• The wavelength centroids are set to the NIST rest 
wavelengths'^, redshifted by the redshift measured 
for Ha. 
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• Each wavelength centroid may vary by up to twice 
the 1 a wavelength uncertainty from the initial fit 
to Ha. 

• All lines must share a common linewidth that may 
vary freely. The initial linewidths are the best-fit 
values obtained from lines in the arc lamps. 

• The flux ratio of the [N 11] dou blet is locked to the 
theore tical value computed by [Storey h, Zeippen 
((2000|. 



Finally, we measure the [S II] emission by fixing the red- 
shift and linewidth to the values obtained for the joint 
Ha and [N II] fit, and fitting the [S II] doublet, allowing 
the flux ratio to vary freely. Table [3] reports the best-fit 
redshifts and flux measurements. Only [S II] detections 
are included. For objects where individual "knots" or 
images were extracted separately, we report the results 
for each component, as well as the summed fluxes over 
all components, and the weighted mean redshift and its 
uncertainty. Figure [2] shows the extracted spectra, the 
1 a error spectra, and the best fit Gaussians to the [N II] 
and Ha emission lines. 

3. PHYSICAL PROPERTIES 

3.1. Lensing mass models 

A crucial component of any study of a lensed galaxy 
is the construction of a mass model for the foreground 
galaxy cluster to derive the lensing magnification and 
demagnify the observed physical properties back to the 
source plane of the galaxy. A detailed description of 
the modeling process for the 8 lens systems considered 
in this paper will be given elsewhere (K. Sharon et 
al. 2012, in preparation). Lensing models for SGAS 
J091538-K382658 and SGAS J134334-f 415509 can also be 



found in Bayliss et al. (2010). To summarize, the lensing 



models are cons tructed using the publicly-available soft- 
ware LENSTOOL ( [Julio et a l.|2007|) with MCMC minimiza- 
tion in the source plane. The clusters are typically mod- 
eled with a Navarro, Frenk & White (1996) profile or a 
Pseudo-Isothermal Ellipsoid Mass Distribution (PIEMD) 
to represent the cluster mass. Additional PIEMDs, 
scaled with luminosity, represent the lensing effect from 
cluster- member galaxies (see Limousin et al. 2007) . The 
redshifts and positions of the multiply-imaged lensed 
background galaxies are used as constraints, and when 
available, the velocity dispersion of the cluster and low- 
uncertainty photometric redshifts of secondary arcs are 
used as priors. 

The source plane reconstruction is done by ray-tracing 
the pixels in the observed frame through the lensing po- 
tential of the best-fit lens model. The ratio of the ob- 
served image size and the size of the model-reconstructed 
source provides an estimate of the magnification. The 
magnification uncertainty is derived from a simulation 
in which many lens models are computed, with parame- 
ters drawn from a set of model parameters from steps in 
the MCMC that correspond to 1 cr uncertainty in the pa- 
rameter space. The magnifications and their associated 
uncertainties for the sample of lensed galaxies are noted 
in Table m 

^ http:/ /www.pa.uky.edu/^peter/atomic/ 



3.2. Stellar mass 

Stellar masses are determined by fitting spectral energy 
distribution (SED) models to the available bro adband 
photometry. Following the recipe presented in |Wuyts| 
et al. (2010), we use t he newest Bruzual and C hariot 
SEDmo3eIs (CB07, see |Bruzual fc Chariot' (^200^)) with 
a Chabrier initial mass funct ioii ( Chabricr 2003 ) and a 
Calzetti dust extinction law (fCalzetti et al. 2000). The 
metallicity is fixed to the value derived from the rest- 
frame optical spectroscopy. To avoid unphysically young, 
luminosity-weighted solutions, the age of the stellar pop 
ulation is required to be > 50 Myr, which rou ghly corre 
spond s to the dynamical timescale at 2 ^ 2 (Erb et al. 



2006c). We have carried out the SED fitting with con 



stant star formation (CSF) models as well as a range of 
exponentially declining star formation histories (SFH), 
parametrized as SFR(t) - e~*/^ with t=10, 50, 100, 
200, 500, 1000, 2000 and 5000 Myr. The data gener- 
ally can not constrain the shape of the SFH; we found 
that most SFHs provide acceptable fits for most sources. 
However, the choice of SFH can have a large influence on 
the stellar population parameters reported by the SED 
fit, especially the d ust extinction and cu rrent star forma- 
tion rate. Recently, Wuyts et al. (2011) have shown that 
when using exponentially declining SI^'Hs, a minimum e- 
folding time r = 3 Gyr is required to match SED derived 
SFRs to SFRs derived from the bolometric UV-I-IR lumi- 
nosity. They find that although shorter e-folding times 
often report lower values, these models will underpre- 
dict the current star formation rate and underestimate 
the age of the stellar population. We similarly find that 
the shortest e-folding times r — 10—50 Myr are often pre- 
ferred in a least-squares sense for our lensed galaxies, but 
the inferred current SFRs are much lower th an th e SFRs 
we derive from the Ha emission flux (see ^3.3). Since 



for the relatively young ages of our lensed galaxies (50- 
300 Myr), there is little difference between a CSF model 
and an exponentially declining SFH with t = 3 Gyr, we 
choose to limit our SED fits to CSF models. We conser- 
vatively include the influence of a declining SFH in the 
systematic uncertainties of the stellar population param- 
eters. Overall, stellar mass is the most robust parameter 
that can be derived from SED fitting, with an average 
fractional uncer tainty of aM /A'h = 0.4 (Shapley et al. 
2005| [Wuyts eTal. 2007 



One important limitation of single component constant 
or declining SFHs is their insensitivity to underlying old 
stellar populations. For a subsample of bright lensed 
galaxies with similar proper ties to the sample s tudied in 
this paper, we have shown in Wuyts et al. ( 2012 1 that the 
contribution of an underlying old stellar population to 
the galaxy's stellar mass is limited to less than twice the 
mass derived from a single component SFH. We expect 
a similar result here and conclude that the stellar mass 
derived from single-component SED fitting provides a 
reasonable estimate of the total stellar mass of the lensed 
galaxies. 

The best-fit values of the stellar mass are reported in 
Table [4j They have been corrected for the lensing magni- 
fication. The statistical uncertainties on the stellar mass 
are derived from 1000 mock realizations of the observed 
SEDs consistent with the photometric uncertainties. We 
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Table 3. Measured redshifts and emission line fluxes. 



source 


rcdshift 


f(HQ) 


f([N II]A6585) 


f([S II]A6718) 


f([S II]A6732) 


SGAS J085137+333114 


1.69320 ± 0.00020 


247 ± 27 


67 ± 10 


36 ±8 


25 ± 31 


SGAS J091538+382658 


1.50200 ± 0.00002 


170 ± 5 


13 ± 10 






SGAS J103842+484926 


2.19678 ±0.00001 


1220 ± 8 


40 ± 11 


68 ±9 


67± 17 


SGAS J103844+484914 












...knotl 


0.96540 ± 0.00010 


63 ± 11 


11 ± 3 






...knot2 


0.96457 ± 0.00006 


120 ±9 


17 ± 5 






...total 


0.96479 ± 0.00005 


183 ± 14 


28 ±6 






SGAS J111124+140901 


2.13855 ±0.00002 


519 ± 7 


34 ± 7 


15 ± 7 


26 ±6 


SGAS J113809+275445 












...knotl 


1.33224 ± 0.00009 


21 ± 2 


1.3 ± 1.3 


3.9 ± 1.8 


8±4 


...knot2 


1.33240 ± 0.00010 


13 ± 2 


-0.9 ± 1.2 






...total 


1.33231 ± 0.00007 


34 ± 3 


0.4 ± 1.8 






SGAS J113809+275439 












...knotl 


0.90905 ± 0.00001 


127 ± 2 


4.4 ± 1.9 


14 ± 2.6 


5± 1 


...knot2 


0.90907 ± 0.00001 


87 ± 2 


2.6 ± 1.7 


8±2 


5± 1 


...total 


0.90906 ± 0.00001 


214 ± 3 


7.0 ±2.5 


22 ± 3 


10 ± 1.4 


SGAS J120924+264052 


1.02057 ± 0.00002 


51 ± 2 


6.7 ± 1.3 






SGAS J134334+415509 


2.09224 ± 0.00005 


155 ± 10 


23 ±4 


19 ±4 


4±6 


SGAS J142039+395454 


2.16256 ±0.00008 


216 ± 13 


34 ± 11 







Note. — Flux units are lO"^'^ erg s~^ cm~^. For objects where we separately extract discrete "knots" or images, we report measurements 
for each component, as well as the summed fluxes over all components, and the weighted mean redshift and the uncertainty in the mean. 
Only [S If] detections are included. 



include a systematic uncertainty related to the choice of 
SFH by repeating the fits to the 1000 mock realizations 
of the observed SED with a range of exponentially de- 
clining SFHs. 

3.3. Star formation rate 

Table [4] reports the current star formation rates for 
our lensed galaxies, corrected for the lensing magnifica- 
tion, as derived from the range of best-fit SED models. 
We compare these to star formation rates estimated from 
the extinction-corrected Ha emission flux. For the lat- 
ter, we have to aperture correct the flux received through 
the narrow NIRSPEC slit to the total intrinsic flux. To 
this end, we smooth the ^-band image of the source with 
a Gaussian to mimick the NIRSPEC seeing conditions. 
The g-band image is used because the rest-frame UV 
emission can be expected to most closely match the Ha 
emission. From the smoothed g-band image, we measure 
the fraction of the total flux (within the 2*FWHM pho- 
tometric aperture) that falls inside the NIRSPEC slit. 
The correction factors are reported in Table [4] The un- 
certainties on this aperture correction are estimated from 
the uncertainty in the seeing (wh ich was quite variable 
during the observations, see ^2.2 1, and a pointing error 
up to 0.25". 



We derive a star formation rate fro m the aperture cor - 
rected Ha flux using the relation from Kennicutt (1998), 
converted to the Chabrier IMF by dividing by a factor 
1.7. We correct for extinction using the Calzetti dust 
extinction law and the reddening E{B — V) derived from 
the best-fit SED models. We assume that the reddening 
of the nebular emission lines is the same as the reddening 
of the stellar light as derived from the SED fit. 



studied in Erb et al. (2006b, 2006c). They use a similar 
strategy of assuming constant star formation histories 
for their SED fits (except when exponentially declining 
models produce a much improved fit - those galaxies are 
shown with open symbols in Figure p|, and dust correct- 
ing the Ha flux with E{B—V) derivea from the SED. Our 
lensed galaxies are broadly consistent with their sample, 
and show a similar scatter of 0.3 dex in the relation be- 
tween SFRua and SFRsed- The open symbols in this 
plot illustrate how SED models with exponentially de- 
clining SFHs with short e-folding times can significantly 
underestimate the current SFR, as discussed in |3.2[ 



The Ha star formation rates are corrected for the lens- 
ing magnification and reported in Table [4] In Figure |3] 
we compare SFRhc to SFRsed- The plot also shows 
the sample of UV-selected star-forming galaxies at z ~ 2 



It is worth mentioning that the lensing magnification 
of our galaxies introduces an additional source of uncer- 
tainty in SFRua- When we derive the total intrinsic Ha 
flux by aperture correcting the flux to the whole galaxy 
and then correcting for the average magnification, we im- 
plicitly assume that the magnification is uniform across 
the galaxy. If the area covered by the slit is more highly 
magnified than the rest of the galaxy, this procedure will 
overestimate the intrinsic Ha flux and vice versa. We 
can avoid this uncertainty by mapping both the galaxy 
and the slit back to the source plane, and measuring 
the aperture correction factor intrinsically in this plane 
(K. Sharon et al. 2012, in preparation). Given the over- 
all agreement with the non-lensed sample from Erb et 
al. (2006b, 2006c) seen in Figure [sj we do not expect spa- 
tial magnification variations to cause a significant bias in 
SFRhu for our lensed galaxies. 

3.4. Electron density 

The doublet ratio of [S II]A6718-A6732 constrains the 
electron density of the nebular gas. We use the IRAF 
task stsdas. analysis. nebular. temden to calculate for 
the four lensed galaxies where the lines are sufficiently 
well detected, assuming an electron temperature of — 
10^ K. The results are reported in Table [s] We com- 
bine these new measurements of electron density with 
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Fig. 2. — Fluxed ID NIRSPEC spectra of the lensed galaxies. The Y-axis shows flux in units of 10"^'^ erg s"'^ cm~^. Spectra are shown 
in black, 1 a error spectra in gray, and the best-fits to the Ha and [N II] emission lines in blue. The wavelength extent of each best-fit 
shows the region over which the continuum level was fit. 



results from the lit erature for three l ensed ga laxies from 
our previous work ( Rigby et al.|2011| ;W uyts et al.|2012 ) 



two lensed galaxies from iHainl ine et al.| ( )2009| ), and one 
et ai.'(2Ul(j|). i'he measurements 



lensed galaxy from Bian et al.'(2Ul(jp 
span the full range of plausibility from the low density 
limit to Tie > 10** cm~'^, which implies a huge range in 
physical conditions in the nebular regions of star-forming 
galaxies at z ^ 2. There is no evidence for the existence 
of an "average" electron density. This issue should be 
investigated further by obtaining density measurements 
for a larger sample of galaxies, as well as deeper spectra 



to more precisely constrain densities for galaxies with ex- 
isting data, since most of the current measurements are 
quite uncertain. 

3.5. Metallicity 

The only direct way to determine the abundance of 
metals in a galaxy's gas reservoir is through a measure- 
ment of the electron temperature Tg. Higher metallic- 
ity results in more cooling via metal emission lines and 
thus lower temperature. However, the auroral emission 
lines (transitions from the second lowest to the lowest 
excited level) required to measure Te, particularly the 
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Table 4. Physical properties of the sample 



source 




log{M,/M0) SFRsED E(B-V) aperture SFRhc 12 + log(0/i?) 






Mq yr~-^ correction Mq yr"-*- 



SGAS 


J085137+333114 


27+1 


9 


c; + 0.2 
"^-0.2 


32+12 





33±0 


14 


2.7±0.4 


(-Q+33 


fl 07+0.08 
»-»'_0.07 


SGAS 


J091538+382658 


25+f 


8 


4+0.8 
^-0.2 


4.0+1 





17±0 


12 


2.0±0.6 


15+23 


o 49+0.28 
^■4^-0.30 


SGAS 


J103842+484926 


9+1 


9 


g+0.3 
="-0.1 


183+_lf 





27±0 


02 


2.9±0.7 


1530+lfo 


O 11+0.11 

8-ii-o.ii 


SGAS 


J103844+484914 


9+2 


9 


Q+0.3 







14±0 


11 


3.2±0.9 


22+_\l 


8-OO-0.08 


SGAS 


J111124+140901 


14+17 


9 


9+0.6 
^-0.4 







05±0 


06 


2.0±0.6 


130+_ll° 


o 07+0. 07 
-0.07 


SGAS 


J113809+275445 


12+2 


9 


r, + 0.3 

^-0.3 


3711 





30±0 


16 


4.1±0.9 


14+1 


,1+0.14 
8-11-0.15 


SGAS 


J113809+275439 


14+2 


8 


7+0.1 
'-0.1 


7+2 
' -2 





17±0 


11 


2.9±0.8 


14+^ 


7.67l°:re 


SGAS 


J120924+264052 




8 


f;+0.3 

''-0.3 


4li 





22±0 


14 


3.5±1.2 




8.60l°:°? 


SGAS 


J134334+415509 


40i* 


9 


4+0.2 
^-0.2 


13+^ 





llitO 


05 


3.5±1.4 


291" 


8.65l°o:g^ 


SGAS 


J142039+395454 


7+2 
' -1 


9 


q+0.2 
^-0.2 


135li 





28±0 


05 


3.7±0.4 


440l-« 


8.671°-? 



Note. — Columns are source name, approximate lensing magnification, intrinsic stellar mass derived from the SED fit, intrinsic SFR 
derived from the SED fit, reddening derived from the SED fit, slit apertur e correction, intrinsic SFR derived from the dust-corrected Ho 
emission, and metallicity derived from the [N II]/Ha ratio as calibrated by|Maiolino et ar]||2008ll 
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Fig. 3. — Comparison between the SFR derived from the SED 
fit and the SFR derived from the Ha emission. The gray data- 
points correspond to the sample of UV-selected galaxies at 2 ^ 2 
studied by Erb et al. (2006b, 2006c). SED models that asssume an 
exponentially declining SFH are distinguished as open circles. 



Table 5. Electron density measurements 



source 


z 


rie (cm 3) 


reference 


SGAS J103842+484926 


2.20 


600+1''°° 


this work 


SGAS J111124+140901 


2.14 


4500l;^50 


this work 


SGAS J113809+275445 


1.33 


160001^750 


this work 


SGAS J113809+275439 


0.91 


low Tie limit 


this work 


RCSGA0327 


1.70 


235+^° 


Rigby2011 


SGAS1527 


2.76 


400+2'^° 


Wuyts2012 


SGAS1226 


2.92 


low Tie limit 


Wuyts2012 


Clone 


2.00 


1270-2540 


Hainline2009 


Cosmic Horseshoe 


2.38 


320-1600 


Hainline2009 


J0900+2234 


2.03 


11 00+^5°°° 
11UU_700 


Bian2010 



most widely used line [O IIIJA4363, become extremely 
weak at metallicities above ~0.5 solar. Beyond z ^ 1, 
the [O IIIJA4363 line has b een detected in only one low- 



metallicity lensed galaxy (Yuan & Kewley 2009) and 



an upper limit has been derived for a highly-magnifie d 
lensed galaxy of average metallicity ( |Rigby et al.pOll ). 
This forces us to rely on more easily detected optical 
emission lines, such as hydrogen recombination lines (Ha, 
H/3) and coUisionally excited forbidden lines ([N IIJA6585, 
[O II]A3727, [O III]A4959,5007). Photoionization models 
are used to calibrate flux ratios of these lines as "strong- 
line" metallicity indicators. Unfortunately, these mod- 
els are subject to significant uncertainties and system- 
atic effects, which result in offsets up to 0.7 dex between 
different strong -line indicators. From a large sample of 
SDSS galaxies, Kewley & Ellison (2008) present recipes 
to eliminate these ott'sets by converting different metal- 
licity indicators to a common calibration scale. 

Using local metallicity calibrations at high redshift as- 
sumes that physical conditions such as ionization param- 
eter and electron density that go into the photoioniza- 
tion models on which these calibrations are based, do 
not evolve with redshift. There is growing evidence sug- 
gesting that this is not a valid assumption. High red- 
shift star-forming galaxies are offset with respect to the 
local star-forming sequence in the diagnostic BPT dia- 
gram of [O in]A5007/H/3 v ersus [N n]A6585/Ha ( [Erb 
et al. 2006a Liu et aL]|2008 ) . Combining measurements 
for six lensed galaxies at 2 ~ 2, we have shown previ- 
ously that their ionization parameter is roughl y twice as 
high as what is measured for local galaxies (Wuyts et 



al.||2012[ ). A robust in-situ calibration of commonly used 
strong-line metallicity indicators at z ~ 2 through direct 
detection of the [O ni]A4363 is crucial to interpret fu- 
ture abundance measurements at high redshift with the 
next generation of telescopes. The increased flux levels 
of lensed galaxies could play an important role here. 

For this analysis, we choose to derive metallicities from 
the flux ratio of [N H] A6585 to Ha, the N2 indicator. The 
close proximity of the lines allows them to be observed 
in the same fllter, thus eliminating flux calibration is- 
sues, and makes their ratio insensitive to dust redden- 
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ing. This is an important point, given that extinction 
is notoriously difficult to estimate robustly at high red- 
shift. Given these benefits, the N2 indicator has been 
used most widely for m etallicity studie s at z > 1. We 
use the calibration from [Maiolino et al.| ( |2008 1 based on 
a joint fit to a sample of local low-metallicity galaxies 
with measurem ents of the [O IIIJA4363 auroral line (Na- 
gao et al.|2006 l and a sample of SDSS DR4 galaxies with 
metal licities determined f rom t he photoionization mod- 
els in 'Kewley & Dopita ( 2002 1 . The calibration has a 
systematic uncertamty of 0.1 dex. It results in metallici- 
tie s ^0.1 dex higher tha n the N2 indicator as calibrated 
by Denicolo et al. ( 2002^ and ^ 0.2 dex higher than the 
calibration by Pettini & Pagel ( 2004 1 . At high metallic- 
ities, the [N IIJ line starts to dommate the cooling and 
saturates, rendering the N2 indicator largely i nsensitive 
to metallicit y at 12 -I- \og{0/H) > 9.1 (in the |Maioliiio 
et al.| ( |2008[ ) calibration). 

The metallicity measurements for the lensed galaxies 
are listed in Table |4] The quoted uncertainties only re- 
flect the statistical uncertainties in the [N II] and Ha 
emission line fluxes. Apart from the systematic uncer- 
tainty introduced by the choice of metallicity calibration, 
one should keep in mind that an integrated metallicity 
measurement from a long-slit spectrum might not be rep- 
resentative of the true galaxy abundance in the presence 
of strong metallicity gradients. 

4. ADDITIONAL DATA FROM THE LITERATURE 

To complement our sample of lensed galaxies at 0.9 < 
z < 2.5, we have searched the literature for galaxies 
within the same redshift range with published measure- 
ments of metallicity, stellar mass, and SFR. We limit our- 
selves to metallicities measured from the N2 indicator. 
Including other metallicity indicators such as the R23 in- 
dex (based on [O II]A3727, [O III]A4959,5007 and H/3) 
would increase the sample, but the uncertainties related 
to the intercalibration of different metallicity indicators 
would severe ly complicate the analys is; the conversions 
presented by Kewley & Ellison ( 2008 1 are derive d lo cally 
and it is not at all clear given the discussion in §3.5| that 
these offsets will be the same at 2; ~ 2. The literature 
sample contains 60 galaxies, inclu ding 6 galaxies with 
upper limits on the [N II] emission dShapley et al.|[2004' 



200F;Torster Schreiber et al.|2006[ 



al., 2009; Qucyrel et al. 



eraIp009F|Bian et al, 
et al. 2011). 



L uet al."2008 ;'Lawet 



2009[ [Wright et al. 2009 ; Q ucyrel 
20T0I [Richard et al.||Blirilligby 



We have converted the values of metallicity, stellar 
mass and star formation rate reported in these studies to 
be consistent with the measurement methods described 
in ^|3j Stellar mass is the most robust parameter that can 
be derived from the SED modeling procedure; different 
assumptions made in terms of the population synthesis 
models and the range of allowed star formation histories 
in the various studies will generally not affect the result- 
ing stellar mass by more than a fac tor of aM,/M^,—OA 
( |Erb et al.[[2006c| [Wuyts et al.|2007[ ). We include this as 
a systematic uncertainty. I'br the star formation rates, 
we include a factor of 2 to aperture correct the published 
Ha flux from the spectroscopic slit to the whole galaxy. 
This is not necess ary for the studies that us e integr al field 
spectroscopy (jForster Schreiber et al. | 2006;|Wright et al. 
20091 [Queyrel et al.[|2009| [Law et al.[[2TO| . The tih'k is 



corrected for dust extinction based on the Calzetti dust 
extinction law and the reddening E{B — V) reported by 
the SED fit. The reddening depends significantly on the 
assumptions made in the SED modeling, even though 
this is not always adequately discussed in the papers. To 
account for this, we assume a systematic un c ertain ty of 



Shapley et aLj ( |2005[ ) and 



(TE(B-v)/ EiB - V)=0.7. In 

Liu et al.| ([2008[ ) the reddening from the SEL^ fiF was not 
reported, but we use their measurement of the Balmcr 
break instead. 

5. RESULTS 

5.1. The relation between stellar mass and SFR 

To situate our sample of lensed galaxies with respect 
to the literature sample, it is informative to place them 
in a diagram of SFR versus stellar mass. Both in the 
local Universe and at high redshift, star-forming galaxies 
follow a tight relation between their rate of star formation 
and their assembled stellar mass, refered to as the main- 



sequence of star formation ( Noeske et al.[[2007 Elbaz et 
ar |2007 Daddi et al.[ 2007). This suggests a continuous 
mode of galaxy growth over extended timescales, p ossibly 
fueled by the accretion of gas through cold flows (Dekel 
eral]|2009 ). 

Figure^ shows that the lensed galaxies (red data- 
points) probe lower stellar masses than have so far been 
studied: the median stellar mass of — 10^'^ Mq is 
a factor of 11 lower than the median stellar mass of the 
literature sample. The current rate of star formation in 
the lensed galaxies is high (median SFR — 29 Mq yr~^), 
which can be explained by the fact that we have selected 
the sample partly based on high rest-frame UV surface 
brightness (see On average, the lensed galaxies lie 
a factor of 10 abov e the main-sequence at z ~ 2 from 
Daddi et al. (20071. This means that their specific star 
formation rate (sdFR — SFR/M^) is an order of mag- 
nitude higher compared to main-sequence star-forming 
galaxies at z ~ 2. 

Strongly star-forming outliers to the main-sequence of 
star formation are typically interpreted as merger-driven 
starbursts (e.g. Rodighiero et al. 2011). Roughly one 
third of star-forming galaxies at z ~ 2 show evidence 
for a recent merger/interaction in spatially resolved ob- 
servations of the gas kin ematics with near-infrared in- 



tegral field spectroscopy ([Forster Schreiber et al 



Law et al.[[2009| [Wright et al.[[2009 



[2009 

I'hese studies have 

mostly targeted the massive end (M* > 10^-^ Mq) of the 
main-sequence at 2; ~ 2 and there is as of yet no evi- 
dence for a correlation between kinematic classification 
and locati on in the stellar mass-SF R plane (see Figure 
16 in Forster Schreiber et al. (2009)). Our lensed galax- 



ies ofier a unique opportunity to extend these kinematic 
studies both to lower stellar masses and more strongly 
star-forming outliers of the main-sequence. 

The lack of galaxies at the low mass - low SFR end 
of the main-sequence is due to selection effects. Both 
our sample of lensed galaxies and the literature sample 
can be described as SFR-limited, since the galaxies are 
mainly selected based on either UV surface brightness 
or Ha emission. As a comparison, th e UV-selec t ed star - 
forming galaxies at z ~ 2 studied by Erb et al. (2006b!), 
shown as small gray datapoints in Figure|4j span a similar 
range in SFR and stellar mass. The full main-sequence 
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Fig. 4. — The relation between SFR and stellar mass. The liter- 
ature sample is color-coded by reference; lensed galaxies from this 
paper are shown in red. The light gray bar indicates the main se- 
quence ( MS) of star fo rmation and its standard deviation as iden- 
tified by p addi et al .l (20071 at z ~ 2. The dashed hne marks 
sSFRs increased by a factor of 10 compared to the main-sequence. 
Small gray datapoints correspond to UV-selected galaxies at 2: ~ 2 
studied by |Erb et al.| | |2006b| . 



only becomes apparent for a near-IR selecte d sample, 
w hich is mass-limited instead of SFR-limited (Daddi et 
al] [2007l ). 



5.2. The relation between stellar mass and gas-phase 
metallicity 

Figure [5] shows the relation between gas-phase metal- 
licity and stellar mass for the combined sample of our 
lensed galaxies and the data compiled from the litera- 
ture. We have divided the sample into two redshift bins 
with roughly equal numbers of sources: 34 galaxies at 
z < 2 (median z = 1.4) and 36 galaxies at z > 2 (median 
z = 2.2). In both redshift bins, the lensed galaxies probe 
lower stellar masses than have been previously studied. 
As a reference, we show th e best-fit mass-metallic ity re- 
lations (MZR) derived by |Maiohno et al.| ( |2008[ ) (con- 

0.07 basM^ SDSS 
0.7 based 
2.2 



verted to the Chabrier IMF 



at z 



data from Tremonti et al. ( 2004[), at z 
on data from Savaglio et al.| (|2005p, and at z 



based o n stacked star-forming galaxies from |Erb et aT] 
(2006a I. The relations are extended as dotted lines to 



stellar masses below the range probed by the respective 
datasets. 

We investigate redshift evolution of the MZR by fit- 
ting linear relations to the data in each bin. Offsets in 
metallicity are more easily quantified with a simple linear 
fit compa red to the polynomial fits used by Mai olino et| 
al. (2008), especially given the limited number of data- 



points and the large scatter. We use the ASURV package 
(revision 1.3; LaValley et al. 1992) to incorporate both 
metallicity measurements and upper limits. The disper- 
sion around a linear fit to the data in each redshift bin is 
0.23 and 0.25 dex re spectively, simil a r to t he dispersion 
of 0.2 dex found by Savaglio et al. (2005) at z ~ 0.7, 
and twice as large as the dispersion of 0.1 dex for SDSS 



galaxies (Tremonti et al. 2004). The gray areas in Fig- 
ure m show the 1 a confidence intervals of the fit when 
we fix the slope to the best-fit slope derived for all data- 
points combined. The y-axis intercepts are 6.10±0.04 at 
z < 2 and 5.94±0.04 at z > 2 respectively, which trans- 
lates into an evolution of 0.16 ± 0.06 dex in the MZR 
between b oth redshift bins . The best-fit polynomial rela- 
tions from Maiolino et al. ( 2008 ) find a similar evolution 
between z ~ 0.7 and z ~ 2.2 at high stellar masses. 

The lensed galaxies studied in this paper provide the 
first direct probe of the MZR at z = 1 — 2 at low stellar 
masses. We derive best-fit linear relations for the sub- 
sample with stellar mass < 10^'^ M^d, which corre 



sponds to the lowest mass bin for which Erb et al. ( 2006a ) 
was able to measure a metallicity for their stacked sam- 
ple at z ~ 2.2. The best linear fit shifts downwards 
by 0.04 ± 0.14 dex in metallicity from the low to the 
high redshift bin (dashed lines in Figure [5]) . Thus, we 
find less redshift evolution of the MZR at low stellar 
masses. A larger sample of metallicity measurements for 
low-mass galaxies is needed to reduce the statistical un- 
certainties, which are large and formally consistent with 
no change in redshift evolution from low to high stel- 
lar masses. However, our measurements definitely do 
not agree with the much larger redshift evolution (up to 
0.4 dex at M^, = 10^ M©) pred icted by extrapolation s of 
the polynomial relations from Maiolino et al. (2008) to 
stellar masses below the measurements at z = 2.2 from 



Erb et al. (2006a 



The mass-metallicity relation can be extended to even 
lower stellar masses when looking at dwarf galaxies. In 
Figur e [6| we show resu lts for 25 nearby dwarf galaxies 



from Lee et al. 



(2006). These authors derived stellar 
masses from the 4.5 /im luminosity and metallicities from 
the [O IIIJA4363 emission line and find a best-fit mass- 
metallicity relation with a scatter of 0.12 de x, similar to 
the scatt er in the relation for SDSS galaxies (Tremonti et 
al.|2004 ) . In the region where the s ample of dwarf ga lax- 



les overlaps with the SDSS sample, Lee et al. (2006) re- 
port a metallicity offset of 0.2-0.3 dex and attribute this 
to an overestimate of the SDSS metallicities at low stellar 
masses due to the limited central region of nearby galax- 
ies probed by the SDSS fibres and/or an offset between 
bright-line metallicity indicators and the [O IIIJA4363 
emission. Figure [6] shows an overall agreement between 
our lensed galaxies at 0.9 < z < 2.5 and the nearby dwarf 
galaxies with comparable stellar masses. However, given 
the different methods used to derive both stellar mass 
and metallicity, there could be significant systematic un- 
certainties. A direct comparison requires measurements 
of the [O IIIJA4363 line at high redshift, which to date 
has been done for on ly on e low-metallicity galaxy (Yuan 
& Kewley 2009, see 01^ . 

5.3. The fundamental metallicity relation 

Given that star formation is responsible for the build 
up of the stellar mass and metal content of galaxies, 
it is instructive to study the dependence of the mass- 
metallicity relation on star formation rate. Several stud- 
ies have reported a correlation between metallicity and 
(s)SFR for SDSS galaxies, where at a fixed stellar mass, 
more actively star-forming galaxies show lower gas-phase 
metallicities (Ellison et al. 2008 Mannucci et al. 2010 
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Fig. 5. — Mass-metallicity relation in two redshift bins. The literature sample is color-coded by re ference; our lensed gala xies are shown 
in red. The solid lines show the best-fit polynomial relations at z = 0.07 , z = 0.7, and z = 2.2 from|Maioli no et al.|l|2008[ l based on data 
from [rVemonti et al. (2004), Savaglio et al. (2005), and Erb ct ah] | |2006a^ ; the dotted lines show their extension to stellar niasses below the 
range probed by the respective datasets. The vertical arrow in the top left corner shows a systematic uncertainty of 0.2 de x, ba sed on the 
fact that the Maiolino et al. (2008) calibration returns metallicities that are ^ 0.2 dex higher than other N2 indicators (see i ]3.5[ ). The gray 
areas correspond to 1 cr confidence intervals of a linear fit to the data with the slope fixed to the best-fit slope derived for all datapoints 
combined. The dashed lines show the best linear fits for the subsample at Af» < lO'^*'^ Mq in each redshift bin. 
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Fig . 6 . — Mass-metallicity relation for our sample of lensed galax- 
ies at 0.9 < 2 < 2.5 (red circles), the sample compiled from the 
literature in the same redshift range (gray circles and upper limits) 
and the sample of local dwarf galaxies from|Lee et al.l|f2006) (black 
squares). The linear fit to the dwarf galaxie s and the best-fit rela- 
tion for the SDSS galaxies from |Maiolino et~al. ( 2008 ) are shown as 
solid black lines. The dashed line corresponds to a Imear extension 
of the SDSS relation to M* = lO^ ^ Mq. 
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ine the dependence ol metaliicity on Sl^'K t 



We exam- 
br our lensed 

galaxies and the sample from the literature in Figure [7] 
For each redshift bin, we divide the combined sample 
into three bins of SFR, chosen to include a roughly equal 
number of sources in each bin, SFR < 20 Mq yr^^, 
20 < SFR < 40 Mo yr-\ and SFR > 40 M© yr-^. 
Figure [7] shows no evidence for a correlation between 



metaliicity and SFR at fixed stellar mass. We quan- 
tify this by computing best-fit linear relations for each 
redshift and SFR bin, keeping the slope fixed as in S|5.2 



For this sample, the highest SFRs do not correspond to 
the lowest metallicities. Given the 1 a uncertainties of 
~0.09 dex on the y-axis intercepts of the fits, their sep- 
aration is not statistically significant; for each redshift 
bin, the best-fit linear relations for the three SFR bins 
are consistent. A recent study of star-forming galaxies 
at z ^ 1.4 with Subaru, limited to high stellar masses of 
> 10®'^ Mq, similarly found only very weak evidence 
for a correla tion between met aliicity and SFR at fixed 
stellar mass dYabe et al.|[20TT] ) 



plies to SDSS galaxies with 



It is worth pointing out t hat the SF R dependence of 
the metaliicity reported by 



Mannucci ct al. (2010) ap- 
ow Si^'Rs ranging from 0.04 



to 6.0 Mq yr~ . The high redshift galaxies studied in this 
paper are much more actively star-forming, with SFRs 
ranging from 1.7 to 1530 Mq yr~^. It would be very in- 
structive, though observationally challenging, to look for 
a correlation between metaliicity and SFR in less active 
high-redshift galaxies with a range of SFRs comparable 
to SDSS galaxies. 



(|2010]) have proposed a fundamental 
R) where galaxy abundance de- 



Mannucci et al. 
metaliicity relation (l-'M 
pends on both SFR and stellar mass, which reduces the 
intrinsic scatter in the local MZR to 0.06 dex. Based on 
median values of the available literature data, they find 
that star-forming galaxies out to z = 2.5 follow the same 
FMR. They explain the observed redshift evolution of the 
mass-metallicity relation through the selection of galax- 
ies with progressively increasing SFRs at hig her redshift. 
In a follow-up paper, Mannucci et al. (20111 extend the 
FMR to lower stellar masses, from > 10^'^ M, 
M, > 



^0 



to 



1O** -^M0, based on --1400 low-mass SDSS galax- 
ie s, ~ 1% of the to t al sam ple of SDSS-DR7 galaxies used 
in Mannucci et al. (20101. The scatter in the relation in- 
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Fig. 7. — Mass-metallicity relation in two redshift bins, color-coded by SFR. The das hed l ines show the best-fit linear relations for each 
SFR bin, keeping the slope fixed to the best-fit slope for all datapoints combined (see j ]5.2| . The vertical arrow in the upper left corner 
corresponds to the systematic uncertainty in the metallicity calibration, as in Figure [s] 



0.5 

To 

■g 
"co 

2 0.0 
DC 



-0.5 



T 




8.5 9.0 9.5 10.010.511.011.5 

log(M*/Mg) 



9.0 



O 8.5 

O 

+ 

- 8.0 



7.5 




• Shapley2004 

• Shapley2005 

• Forster2006 
Liu2008 
Bian2010 
Queyrel2009 
Law2009 
Wright2009 
Richard2011 

• Rigby201 1 

• this work 



0.36 



0.10-? 



0.03 



0.01 



-12 -11 -10 -9 -8 
log(sSFR) 



-7 



Fig. 8.— Residual from the |Mannucci et al.|f2010t FMR. The 
literature sample is color-coded by reterence as m Figure |4] and 
Figurejs] our lensed galaxies are shown in red. The gray area snows 
the 1 a uncertainty in the FMR for SDSS galaxies (Mannucci et 
al. 2010, 2011) 



Fig. 9. — Metallici ty as a function of speci fic star formation rate, 
based on Figure 2 in [Mannucci et ah | |2011| . The gray areas show 
68% and 90% confidence mtervals tor sUbb galaxies. The literature 
sample is color-coded by reference as in Figure |4] Figure [s] and 
Figure [s] our lensed galaxies are shown in red. 



creases dramatically, up to 0.35 dex at = 10 M0. 
Figure |8] shows the residual from the FMR for our sam- 
ple of lensed galaxies and the sample from the literature. 
There is a general agreement at all stellar masses, the 
median residual is 0.29±0.28 dex. The gray area in Fig- 
ure [8] corresponds to the 1 a uncertainty in the FMR 
(Mannucci et al. 2010, 2011) - the high redshift sample 
generally shows a larger scatter. 



Figure |9] is based on Figure 2 in Mannucci et al. (2011 1 
and shows metallicity as a function of spechc star for- 
mation rate. The gray areas contain 68% and 90% of 
the SDSS galaxies used for their derivation of the FMR. 
The scatter in the metallicity increases dramatically with 



sSFR. The high-redshift galaxies studied in this paper are 
generally consistent with the SDSS confidence interval. 
The sSFRs of the lensed galaxies are on average an order 
of magnitude of higher than the literature sample, which 
generally follows the main-sequence of star formation at 
2 (see also § 5.1). The lensed galaxies thus probe a 
1 sSFR, not present in the local SDSS 
metallicity continues to decrease 



new range of hig 
sample. In this range 
with sSFR, and the scatter remains substantial 



6. DISCUSSION 

As discussed in the introduction to this paper, the cor- 
relation between stellar mass and metallicity is a natural 
consequence of the conversion of gas into stars within 
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galaxies, mitigated by gas exchanges with the environ- 
ment through inflows and outflows. Since gas flows are 
hard to observe directly at high redshifts, the observed 
mass-metallicity relation can indirectly constrain their 
relative importance for the chemical evolution of galax- 
ies. In this discussion, we limit ourselves to simple ana- 
lytical models of chemical evolution. This undoubtedly 
does not capture the full complexity and interplay of the 
physical processes involved, but we will show that the 
combination of large scatter and significant systematic 
uncertainties in the key observables, do not allow con- 
straints on even these simple models. 

6.1. Analytical models 

In the closed-box model of chemical evolution, galaxies 
do not interact with their environment and the metal- 
licity is fully determined by the galaxy's gas fraction 
fgas = Mgas/{Mgas + M*) and by the integrated stel- 
lar yield y, the mass of metals produced by massive stars 
normalized to the stellar mass locked-up in long-lived 
stars and stellar remnants. 



Z = VHfgas) 



(1) 



In this equation, Z is the metal mass fraction, the frac- 
tion by mass of elements heavier than helium. The in- 
tegrated yield y is derived from the yields of individ- 
ual stars, weighted with the adopted initial mass func- 
tion. Due to remaining uncertainties in stellar evolu- 
tion and (explosive) nucleosynthesis calculations, the in- 
tegrated yield is not known to better than a factor of 2 
The solid line in Figure [TO] shows the closed-box model 
for y = 0.015 ± 0.005, the appropriat e inte grated yield 



for the Chabrier IMF (Finlator fc Dave 2008 Peeples 



Sz Shankar 2011| ). Simple analytical extensions of the 
closed- box model which incorpora te inflows an d outflows 
are compreh ensiv ely discussed by l Erb ( 2008[ ) - see also 
Pagel| ( |1997[ ) and |Spitoni et al.| ( |20TUf : Dut nows are as- 
sumed to occur at a constant traction of the SFR (fo), 
with the outflowing gas showing the same metallicity as 
the remaining gas reservoir. In this model, the outflow 
rate does not depend on the mass of the galaxy. The 
dashed line in Figure 10 corresponds to a model with 
fo = 2.0, showing a shallower increase of the metallic- 
ity with decreasing gas fraction due to the removal of 
metals by the outflows. To add inflows, the model as- 
sumes that gas is accreted at a constant fraction of the 
SFR (fi) and that the metallicity of the inflowing gas is 
negligible with respect to the metallicity of the galaxy's 
gas reservoir. The inflow rates are generally required to 
be sufficient to roughly replenish the processed gas; an 
accretion rate below 85% of the gas processing rate fails 
to reproduce the observed trends of s tar formation rate 
with galaxy age at z ~ 2 (Erb (2008). The dot-dashed 
curve in Figure [lO] shows the added effect of inflows at a 
rate fi = fo + 0.95. The metallicity is seen to asymptote 
to its equilibrium value Zeq = y/ fi at low gas fractions 
where a balance between the rate of enrichment from star 
formation and the rate of dilution from infall of metal- 
poor gas sets in. Figure [T0| shows that the uncertainty in 
the integrated yield creates significant overlap between 
the different analytical models in the Z-fgas plane. 

6.2. Gas fractions 



A crucial component of the analytical models pre- 
sented above is the mass of a galaxy's gas reservoir. 
(Sub)millimeter observations of molecular gas have only 
recently become feasible for main-sequence star-forming 
galaxies at z = 1 — 2 through significant advances in the 
sensitivity of the IRAM Plateau de Bure millimeter Inter- 
ferometer (PdBI; Guilloteau et al. 1992). Additionally, 
the Atacama Large Millimeter Array (ALMA) will pro- 
vide unique opportunities for enhanced sensitivity and 
spatial resolution for molecular gas studies in the near 
future. Current samples of high-redshift main-sequence 
star-forming galaxies observed with the PdBI are still 
small (15 sources at z = 1 — 1.5 and 10 sources at z ~ 2.3) 
and biased towards the massive end of the main-seq uence 
with M« ^ 3 X 10^° M(D &nd SFR > 40 Mq yr"! ( |Tac- 
coni et al.|| 20To{ |Daddi et al.||2010[ ). 



In the absence ot direct measurements of the gas mass, 
it has become common practice to use the empirical 
Kennicutt-Schmidt (KS) law to estimate the gas surface 
density from the SFR density: ^sfr ^ ^^as (Kenni- 
The half-light radius R is typically used 



cutt 



1998). 



to convert between SFR and SFR surface density, such 

that Mgas = ^S'i^i?l/"i?(2n-2)/n^ rpj^^ ^^^^^ ^ 

calibrated from the current small samples of gas mea- 
surements for main-sequence star-formin g galaxies at 
z = 1 — 2; results vary from n = 1.17 (Genzel et al. 



2010) 



2010| ) ton = 1.42 paddi et al. 

estimate the gas masses for our sample of lensed 
galaxies, we compute the half light radius R by map- 
ping the photometric aperture that contains half of the 
total flux back to the source plane of the galaxies using 
the lensing model. For the sample from the literature, 
we use the published half-light radii when available - no 
measurements were publish ed for the g alaxies taken from 



Shapley et al. ( 2005) JLiu et aL] pOOSl ) , |Bian et al.J pOlO I 
and |Rigb°et al.^ ( ,201l) , which c orresponds to 42% of the 
sam]Dle. We use the KS law from Daddi et al. ( 2010 ) with 



1.42 and^ = 4.98 x 10*^; t 



Daddi et aL ( 2010[ ) with 
re best-fit KS law from 
Genzel et al. (|2010| with n — 1.17 leads to gas fractions 
that are on average 15% larger. 

6.3. Observational constraints on analytical models? 

Figure [T0| shows the observed metal mass fractions Z, 
and the gas fractions fgas infered from the KS law for our 
lensed galaxies and the subsample from the literature for 
which the half-light radius has been measured. The data 
span almost the full range of gas fractions and there is no 
obvious correlation with the metal mass fraction. The- 
oretically, the phase space to the left of the closed-box 
model is not allowed. At a given gas fraction, a galaxy 
can not have a higher metal mass fraction than predicted 
by the closed-box model, since interactions with the en- 
vironment through metal-rich outfiows and metal-poor 
infiows will only lower the metallicity. An integrated 
yield y = 0.05 is required to obtain a closed-box model 
that serves as an upper limit to the metallicity for mos t 
datapoints (magenta line in the left panel of Figure 10 ) 



Even though the integrated yield is not known to great 
precision, a value oi y = 0.05 is unphysically high. 

The high metal mass fra ctions i n Figure [T0| could al- 
ternatively suggest that the Maiolino et al. ( 2008 ) metal- 
licity calibrat ion is biased high. Wh en applied to SDSS 
galaxies, the Maiolino et al. (20081 calibration finds a 
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Maiolino et al. (2008) Pettini & Pagel (2004) 




Fig. 10. — Relation between metal mass fraction Z and gas mass fraction /gaa. We show a closed-box model (solid line), an outflow 
model with /□ = 2.0 (dashed line), and a model with outflows at fo = 2.0 and inflows at fi = fo + 0.95 (dot-dashed line). In each case, the 
lines correspond to y = 0.015 and the gray areas show the confidence intervals that correspond to a range of yiel ds v = 0.010 — 0.020. The 
horizontal arrow shows the systematic uncertainty in fgas from the uncertainty in the slope of the KS law (see i ]6.2[ l. In the left panel, an 
integrated yield y = .05 is required to obtain a closed-box model that serves as an upper limit to the m etal mass traction as calibrated by 
[Maiolino et al.| | |2008| (magenta line). The right panel shows the calibration by [Pettini fc Pagel| | |200"4] |. 



MZR very similar to the result from ' Tremonti et al.| 
(|2004|). This MZR lies towards the high-metallicity end 



of the . 7 dex spread in metallicity found by Kewley ^ 
Ellison (2008) for 10 different calibrations of the local 
MZR (see their Figure 2). In the right panel of Fig- 
ure |10[ we show metal mass fractions derived from the 
third order polynomia l calibr ation of the [N II] /Ha ratio 



from Pettini & Pagel (2004 1. This returns metallicities 
tha t are on average (J.2 d ex lo wer than the [Maiolino et| 
al. ( 2008 1 calibration (see jj3.5| and falls in the mid dle of 
tEe range of MZRs from [Kewley k EllisonI ( |2008[ ). For 
this calibration, an integrated yield of y = 0.1)15 ± 0.005 
is no longer at odds with most of the data. 

Our lensed galaxies have very high gas fractions of 
70-95% and consistently show larger metal mass frac- 
tions than allowed by th e closed-box model, even for the 



Pettini fc Pagel (|2004) metallicity calibration. Daddi 



et al.[ ( |2010p and |Genzel et al 



(2010) have shown that 



both local and high-redshilt submm galaxies, whose ex- 
treme SFRs are powered by merger-driven starbursts, 
have much higher SFR densities at a fixed gas surface 
density; their KS law has a similar slope compared to 
main-sequence star-forming galaxies, but is normalized 
higher by ^ 0.9 dex. As seen in § 5.1 the lensed galax- 
ies are outliers of the main-sequence of star formation at 
z ~ 2 with 10 times higher specific star formation rates, 
which is typically attributed to merger-driven star for- 
mation. If these galaxies are indeed powered by a recent 
merger, we are overestimating their gas masses by 20% 
on average when using the lower normalization of the 
KS law for main-sequence star-forming galaxies. This 
uncertainty can be addressed with either integral field 
spectroscopy to establish the kinematics of the lensed 
galaxies or direct measurements of their gas mass. 



In summary, the current situation does not allow ob- 



servational constraints on simple anal ytical models of 
chemical evolution introduced in |6.1[ We have iden- 
tified three sources of systematic uncertainty that need 
to be addressed. First of all, a robust absolute calibra- 
tion of the metallicity scale is required. We need detec- 
tions of the [O IIIJA4363 emission line at high redshift to 
provide an in-situ direct metallicity calibration via the 
electron temperature. Until the next generation of tele- 
scopes (e.g., JWST) comes online, the increased flux lev- 
els of lensed galaxies may provide the only path forward 
in this respect. Secondly, we need direct measurements 
of the gas mass for large samples of star-forming galaxies 
at z = 1 — 2, both for main-sequence galaxies and out- 
liers. This is essential to reduce the uncertainty in both 
the slope and the normalization of the KS law. Even 
better, a significant galaxy sample with both metallic- 
ity and gas mass measurements can avoid the KS law 
entirely. Lastly, the uncertainty on the integrated yield 
causes large overlaps between analytical models in the 
Z-fgas plane, limiting our ability to distinguish between 
them. 

7. SUMMARY 

This paper presents measurements of the stellar mass, 
star formation rate and metallicity for 10 lensed galaxies 
at 0.9 < z < 2.5, identified in the SDSS DR7. We sup- 
plement this sample with 60 mostly non-lensed galaxies 
from the literature in the same redshift range, limited 
to sources with metallicities derived from the [N II] /Ha 
flux ratio. Our conclusions are summarized as follows: 

1. The lensed galaxies probe stellar masses that are on 
average a factor of 11 lower than have previously 
been studied in the literature. Their current rate 
of star formation places them on average a factor 
of 10 above the main-sequence of star-formation at 
z ~ 2, which can be explained by the fact that they 
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were selected for near-IR spectroscopy from a much 
larger sample of lensed galaxies (M. D. Gladders et 
al. 2012, in preparation) partly based on high rest- 
frame UV surface brightness. 

2. We find a redshift evolution of 0.16 ± 0.06 dex in 
the mass-metallicity relation (MZR) from z ~ 1.4 
to z ^ 2.2. The poly nomial best-fit MZRs from 



Maiolino et al. 



( 2008) ba sed on data a.t z — 0.7 

from Savaglio et aT "P005| and stacked galaxies at 
z = 2.2 from Erb et ai. ( |20U6a I show a similar evo- 
lution at high stellar masses. The lensed galaxies 
provide the first direct constraints on the MZR at 
z = 1 — 2 at stellar masses below 10^'® M©. Con- 
trary to_tlie_extra£ola best-fit MZRs 
from Maiolino et al. (20081, we find less redshift 
evolution at these low stellar masses. 



3. We do not see a correlation between metallicity and 
SFR for our combined sample, which probes higher 
SFRs than the SDSS galaxies for which this corre- 
lation has previously been reported. We find a gen- 
eral agreement with the local fundamental relation 
between metallicity, stellar mass and star forma- 
tion rate (Mannucci et al. 2010, 2011) at all stellar 
masses, albeit with larger scatter. 

4. We use the Kennicutt-Schmidt law to infer gas 
masses for the sample and investigate whether we 
can constrain simple analytical models of chemical 
evolution with the observed metallicities and gas 
fractions. An integrated yield y — 0.05 is required 
such that the metal mass fractions are bounded by 
the closed-box mode l This is unphysicall y high, 
and suggests that the Maiolino et al. ( 2008 ) metal- 
licity calibration is biased high. 



This paper has shown that strongly lensed galaxies of- 
fer unique opportunities to extend studies of the mass- 
metallicity relation to the low-mass end of the galaxy 
population at z = 1 — 2. We are currently obtaining 



spectroscopy for a larger sample of lensed galaxies to 
reduce the statistical uncertainties and better constrain 
the redshift evolution of the relation at low stellar masses. 
Additionally, we expect the increased flux levels of lensed 
galaxies to play a crucial role in the path towards reduc- 
ing the current systematic uncertainties in the metallicity 
calibration and KS law. 
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